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Abstract 

Because of the infrared renormalons, it is difficult to get power accuracy in 
the traditional approach to the Wilson's operator product expansion. Based 
on a new perturbative renormalization scheme for power-divergent operators, 
I propose a practical version of the OPE that allows to calculate power cor- 
rections to desired accuracy. The method is applied to the expansion of the 
vector-current correlation function in QCD vacuum, in which field theoretical 
status of the gluon condensate is discussed. 
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Wilson's operator product expansion (OPE) |J has become one of the most powerful tools 
in modern field theory. Indeed, many graduate text books on field theory devote an entire 
chapter treating the subject 0]. Direct applications of the OPE in Quantum Chromody- 
namics (QCD) include deep-inelastic scattering, heavy-quark expansion, weak-nonleptonic 
decays, and QCD sum rule calculations, etc. Generally speaking, the OPE provides the 
starting-point for any perturbative calculations in hard processes. The underlying principle 
of the OPE — scale separation and factorization — has spawned many uses of the so-called 
effective field theories in nuclear and particle physics. 

Given the success of the OPE, it may be surprising to realize that the standard text- 
book approach to the expansion is too formal to be used in some applications. A problem 
occurs when one seriously considers the so-called power corrections to the leading term in 
the expansion. To explain, let me first remind the reader some basics of Wilson's expansion 
@. Consider a time-ordered product of two currents separated by, say, an Euclidean short 
distance £ 2 . According to Wilson, the product can be expanded as, 



where, for simplicity, I have neglected possible Lorentz indices. The short distance physics 
above the factorization scale fj 2 is included in the coefficient functions Cj, ordered according 
to the descending singularities as £ 2 — > 0. Oj(/i 2 ) are renormalized operators with increasing 
dimensions, designed to account for long distance physics below /i 2 . The /i 2 dependences 
are cancelled between Cj and Oi term by term in the expansion. 

The expansion in principle can serve to define the composite operators Oi if one knows 
how to classify the short distance singularities of an operator product. In practice, however, 
a reverse procedure is followed in constructing the expansion: A tower of operators with 
right quantum numbers are first chosen, and then the coefficient functions are calculated 
by sandwiching the expansion in a set of perturbative states. Implicitly assumed in this 
standard procedure, though not essential to the principle of the OPE, is that the power 
divergences in the operators Oj are subtracted with a perturbative normal ordering. For the 
operators with vacuum quantum numbers, such as F 2 = F af3 F a/ 3, it means the subtraction 
of the perturbative-vacuum expectation value. For others, it means use of the dimensional 
regularization in which the power divergent integrals / d d k/k a = are taken to be zero. 

It was first pointed out by t' Hooft that the coefficient functions thus obtained, 
Ci{a s ) = J2n a ™ c in, m &y contain the infrared (IR) renormalons — a jargon for c in ~ n\ to 
grow factorially with a fixed sign. Although this has never been proved conclusively in QCD 
except in certain limiting cases such as the number of flavors Nf — > oo, many have taken the 
observation seriously. If true, Ci is a non-Borel summable series and is genuinely infinite. 
If treating Cj as an asymptotic series, it is easy to see the uncertainty in regularizing the 
series, such as truncating it around the minimal term, or regulating the singularity in the 
Borel plane after a Borel transformation, spoils a clear-cut classification of Cj(£ 2 ) in their 
singularities. At practical level, the power corrections to the OPE become ill-defined due to 
the uncertainty in the leading-order coefficient function. 

For applications where power corrections are small, like heavy-quark expansion for the 
top flavor, the problem is of only formal importance. However, there are applications where 
the understanding of power corrections is absolutely essential. In extracting the strong cou- 
pling constant a s from the deep-inelastic sum rules, such as Bjorken and Gross-Llewellyn 
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Smith sum rules, the so-called higher-twist contributions are an important source of theo- 
retical uncertainty ||. In the charm and bottom systems, I/tuq corrections to Isgure-Wise 
symmetry relations [|5] valid in the heavy quark limit are substantial. Finally, the entire 
QCD sum rule phenomenology initiated by Shifman, Vainshtein, and Zakharov || relies on 
a meaningful understanding of power corrections at fundamental level. 

The IR renormalons in QCD are generated from perturbative calculations in infrared 
regions of Feynman diagrams. An indication for this is that the high-order diagrams effec- 
tively renormalize the coupling constant appearing in the lower-orders. The infrared Landau 
singularity in the resulting running coupling is essentially the source for the n\ behaviors at 
high orders J7|. In light of this observation, the appearance of the IR renormalons in the 
OPE is somewhat surprising, for the OPE is devised precisely to account for soft regions 
of perturbative diagrams with the matrix elements of composite operators. The puzzle was 
solved first by F. David ||, who showed that the power- divergent composite operators de- 
fined with the standard subtraction scheme also contain the IR renormalons. And the IR 
renormalon singularities in the coefficient functions and operators matrix elements cancel. 

Although the IR renormalons do not invalidate the standard form of the OPE, they 
pose practical problems in carrying out the expansion to power accuracy. The coefficient 
functions and the power-divergent composite operators are both ill-defined, and one must 
find a consistent way to regularize them. Grunberg |J and Mueller have both suggested 
ways to regularize the perturbation series in coefficient functions. However, computations of 
higher- twist contributions in their schemes appear difficult. On the other hand, in a recent 
work by Martinelli and Sachrajda [|ll|], a non-perturbative renormalization of composite 
operators is proposed. It remains to see, however, how practical the method can be used to 
compute subtraction in coefficient functions. 

Since the IR renormalons arise from soft regions in Feynman diagrams, a physical solution 
to the problem is to avoid perturbative calculations in these regions. The idea is not new, 
as it has been advocated by Novikov et al. since many years ago. However, the only 
serious study in this direction was the work by Mueller ||13|| , which has largely been ignored 
in the literature. The main problems with this approach, to my opinion, have been two: 
first, how to make systematic soft subtraction in the coefficient functions? second, how to 
calculate the non-perturbative matrix elements of composite operators consistent with the 
subtraction. In this paper, I hope to provide answers to both questions through the example 
of the vector-current correlation function in QCD vacuum. 

To begin, let me recall once again the basic principle of the OPE: the composite operators 
are introduced to account for soft contributions in Feynman diagrams. Thus if the OPE 
is applied in a perturbative state, the composite operators shall contain perturbative soft 
contributions. In the standard regularization for the composite operators, however, the 
perturbative contributions are completely subtracted along with the power divergences. This 
is convenient for calculating the coefficient functions, but the operators defined in this way 
cannot account for the perturbative soft contributions in Feynman diagrams for correlation 
functions. Therefore, a fundamental solution to the IR renormalon problems is to find a 
practical way to regularize power divergences of composite operators while maintaining their 
soft contributions to a perturbative calculation. 

Composite operators can be defined entirely through their insertions into Green's func- 
tions |TJ]]. As such, the renormalization of the Green's functions determines the renormal- 
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ization of the composite operators. According to the standard renormalization theory, we 
need to consider only the Green's functions with non-negative superficial degree of diver- 
gence, which is defined for a graph G as, 5(G) = 4 — do — d^n^ , where d Q and d^ are 
the dimensions of the operator O and field <f>, and is the number of external lines in 
the graph. The counter term operators can be constructed from the overall subtraction of 
the primitively divergent Green's functions. Now it is well-known that the renormalization 
subtraction have large freedom with finite contributions (renormalization group). In QCD, 
since the low momentum regions of a Feynman diagram is not perturbatively calculable, one 
shall not subtract divergences along with any soft contributions. In the following, I call such 
subtraction as "minimal" subtraction of power divergences. 

To illustrate a systematic way of doing minimal subtraction, let me construct the com- 
posite operator [0 2 ] in the 4 theory. For our purpose, we are interested in only one insertion 
of the operator. According to power counting, two Green's functions have superficial diver- 
gences, (O|T0 2 0(xi)0(a;2)|O) and (O|T0 2 |O), with superficial degrees of divergence and 2, 
respectively. Therefore, the renormalized composite operators can be defined as, 

[0 2 ] = Z 2 <f) 2 + Z I (2) 

where Zi are renormalization factors. Z 2 is determined by subtraction of primitive log- 
divergent diagrams in (0\T<p 2 <j)(xi)<j)(x 2 ) |0). On the other hand, Z is determined by primitive 
quadratically-divergent diagrams in (O|T0 2 |O). In the standard approach, one writes, Z = 
— Z 2 (O|T0 2 |O) P ert- The operator [<f> 2 ] defined in this way depends on perturbative calculations 
in soft regions. On the other hand, The minimal subtraction can be defined as follows. 
Consider any perturbative diagram for (O|TZ 2 2 |O). Subtract only the contributions from 
regions where all the momenta in the diagrams are larger than some scale A 2 . For instance, 
suppose we have the following contribution to (0\TZ 2 (f> 2 \0), 

roc roc roc 

/ d% d 4 k 2 ... d 4 k n f(h,...,k n ) (3) 
Jo Jo Jo 

Then we subtract the part when all the integrations run from A 2 to oo, 

roc roc roc 

Z = - ... d%d%...d 4 k n f(h,...,k n ) (4) 

J A 2 J A 2 J A 2 

Thus, Zq contains only hard contributions. It remains to see however, the remaining part 
of the Green's functions, with mixtures of soft and hard momentum flows, are finite. 

The proof for the finite remaining part is done by explicit construction. It is easy to see 
through working with two or higher-loop diagrams that 

(O|[0 2 ]|O) pcrt =E(-ir +1 / d%... d'k n 

x(0\TZ 2 <p 2 <p(k 1 )... ( p(k n )^-k 1 - k 2 ... - k n )\0) tr ™ 
x<0|T<Kfci)...<Kfc n M-A;i... - k n )\0)(fi 2 ) (5) 

where the first factor in the integrand (0\TZ 2 <fi 2 <fi(ki)...<fi(k n )(f)(—k 1 ... — fc„)|0) trun is a con- 
nected, truncated Green's function with external momenta k±, k 2 , ... ,k n , —k\... — k n . It is 
finite because all logarithmic sub-divergences are cancelled by the renormalization factor Z 2 . 
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The second factor {0\T<f>(ki)...<f>(k n )<f>(— k\... — fe n )|0) is an ordinary Green's function with 
all internal lines restricted to k 2 > /i 2 . 

To illustrate the use of the minimal-subtracted operators in the OPE, I consider the 
vacuum correlation function of the conserved vector currents in the chiral limit (m q = 
0) ife^d^TlJ^J^O)] = (-g^q 2 + g M ^)n(g 2 ) . Define Adler's ^-function through, 
D{q 2 ) = -{2it) 2 q 2 {<m/dq 2 ). The OPE for D{q 2 ) at the space-like q 2 reads (Q 2 = -q 2 ), 

D(Q 2 ) = C (a s , A 2 ))/ + C,(a s ) 2 ™fg} A2) + ... , (6) 

where J is a unit operator and the renormalization scale /i 2 is chosen to be Q 2 . F 2 (A 2 ) is 
minimally subtracted operator without the quartic divergence, 

F 2 (A 2 ) = F 2 - (0|F 2 |0) pcrt (all k 2 > A 2 ) . (7) 

The renormalized is free of logarithmic divergences. In dimensional regularization and 



with covariant gauge fixing, we have [15 



^ , <91nZ ax ~ , d\nZ„ ,d\nZ 3 1 d\nZ 3 ^ 

F « = (1 + g ^f )F + " A ^r + 2 9 ^r )0gf 

where A is a gauge parameter, and Z g , Z 2 , Z 3 and Z are renormalization constants for the 
gauge coupling g, quark fields ip, gauge potential and ghost field u, respectively. All 
operators on the right-hand-side are unrenormalized and, Oqf = — ^r(c^A-») 2 ~ §ix~A^ 
where S is the gauge-fixed QCD action. 

Let us see how the coefficient function Cq and the non-perturbative matrix element 
F 2 (A 2 ) can be consistently calculated. First, sandwiching the expansion in Eq. (|6]) in the 
perturbative vacuum state, we have, 

C {a s , A ) = C (a s ,(Jj - C 4 (a s J (9) 

where Co(a s , 0) is the coefficient function calculated when assuming the perturbative matrix 
element of F 2 (A 2 ) vanish. The contributions to Cq{oc s , 0) from one or two soft-gluon or ghost 
lines are subtracted by the second term in Eq. (|9[) which, according to our definition, can 
be computed through a formula like in Eq. (|5]). 

Since perturbative diagrams are most conveniently evaluated in the dimensional regular- 
ization and modified minimal subtraction scheme, we shall study the OPE in this scheme 
(labelled by MS). Co(<y s , 0) has been calculated to three-loops JlTj (assuming three flavors), 

u K „ „ ( a s \ 2 „ . f a s ^ 3 



C (a a> 0) = l + -i + 1.6f-^) +6-4(-M + ... , (10) 

and C 4 has been calculated at one- loop 0, C 4 (a s ) = 1 + Zjj* + .... To illustrate the soft 
subtraction, I have calculated the perturbative matrix element of F 2 (A 2 ) up to two-loops, 
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(0|F 2 (A 2 )|0)^ t = ^ (l + (2-41 + 1.25 In (11) 

Finally, we have coefficient function for the unit operator in Eq. (6), 

C™(a s , A 2 ) = 1 + ^ (l - 2^ + 1.67(^) 2 (l - ^ (2.89 + 1.50 In £f) + ... (12) 

The above series is free of the IR renormalon at b = 8tt//3o on the Borel plane. 

To obtain the non-perturbative matrix element of F 2 (A 2 ) (gluon condensate) in dimen- 
sional regularization, we first write, 

(0|F 2 (A 2 )|0) = (0|F 2 |0) - (0|F 2 |0) pcrt + (0|F 2 (A 2 )|0) pert . (13) 



According to the Joglekar-Lee theorems |]16|| , the matrix elements of equations-of-motion 
operators and BRST-exact operators vanish. Thus, the first two matrix elements on the 
right-hand side of Eq. ( |T3| ) depend on the bare operator F 2 only, and they transform 
homogeneously under a scale transformation. Furthermore, the study on the trace anomaly 
of the QCD energy-momentum tensor indicates (3{g)F 2 is renormalization scheme and scale 
independent gS^. Thus, (0|F 2 10>^- (0|F 2 |0)pf t = Z [(0|F 2 |0) LAT - (0|F 2 |0)^] w here 
Z = P LAT j p MS depends on lattice and MS /5-functions. Factorizing Z in Eq. (|I3"D, we find 

<0|F 2 (A 2 )|0p = Z [(0|F 2 |0) LAT - (0|F 2 |0)^ + ^~ 1 (0|F 2 (A 2 )|0) P ^f t J (14) 

The above is a practical definition of the gluon condensate that are renormalon-free and 
consistent with the MS evaluation of the coefficient functions. 

The first two terms in Eq. (O ) are what has been used as a lattice definition of the 



gluon condensate in the literature f2l| . However, (0|F 2 |0)p A ^ is an ill-defined perturbation 
series due to the IR renormalon at b = 87r/f3 . The perturbation series has been numerically 
evaluated up to eight-loops by Di Renzo et al. in the quenched approximation |p2| (see also 
|2j, 



(0|F 2 |0)^T = -Al + 4.01^ + 64.08(^M) 2 + 1321.8 6 (^) 3 + ..) (15) 

To extract the gluon condensate with a good accuracy, one is faced with two opposite 
requirements for the lattice spacing a. To ensure that the calculations are at the continuum 
limit, one shall take as small a as possible. On the other hand, for small a, the cancellation 
of quartic divergences requires extremely good accuracy in the perturbation series. Even 
for a = 0.17 fm (/3 = 5.7), one needs to evaluate the perturbation series to at least one 
percent. High-accuracy is not possible without regularizing the IR renormalon in the series. 
[Although for expansions with the bare lattice coupling, the smallest term is postponed to 
higher orders, the intial decrease of the terms is very slow due to large tadpole contributions.] 
The last term in Eq. (|I4D is introduced precisely for the IR renormalon regularization. Using 
the one-loop relation between the lattice and MS couplings p4j, I find, 



_1 (0|F 2 (A 2 )|0)p^rt = ^ [l + (20.87+ 1.25 In |J + 2.75 ln(QV)) ^ + ... ) (16) 
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To determined the subtraction up to three-loops, we need two-loop relations between lattice 
and MS couplings and the matrix element (0|F 2 (A 2 ) |0)^t up to three-loops, the former 
has recently been calculated by Luscher and P. Weisz [25|| . We wish to make a three- loop 
subtraction calculation in the near future. 

To summarize, I have introduced a minimal subtraction scheme to regularize power 
divergences of composite operators. When these operators are used in Wilson's OPE, the 
coefficient functions are automatically free of the IR renormalon singularities. The matrix 
elements of the composite operators can be calculated in non-perturbative methods such as 
lattice QCD. Thus, within the scheme one can meaningfully discuss power corrections to 
desired accuracy. The scheme can be straightforwardly applied to deep-inelastic sum rules 



where the mixing of higher-twist operators with lower-twist ones is an important issue |26 
It can also be used in the heavy quark effective theory to define the pole masse for the 
heavy-quark expansion | 27fl . This will be discussed in a separate publication. 
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